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Abstract 

We give a representation of the parity-even part of the planar two-loop six-gluon MHV amplitude 
of M = 4 super- Yang-Mills theory, in terms of loop-momentum integrals with simple dual confor- 
mal properties. We evaluate the integrals numerically in order to test directly the ABDK/BDS 
all-loop ansatz for planar MHV amplitudes. We find that the ansatz requires an additive remain- 
der function, in accord with previous indications from strong-coupling and Regge limits. The 
planar six-gluon amplitude can also be compared with the hexagonal Wilson loop computed by 
Drummond, Henn, Korchemsky and Sokatchev in arXiv: 0803. 1466 [hep-th]. After accounting for 
differing singularities and other constants independent of the kinematics, we find that the Wil- 
son loop and MHV-amplitude remainders are identical, to within our numerical precision. This 
result provides non-trivial confirmation of a proposed n-point equivalence between Wilson loops 
and planar MHV amplitudes, and suggests that an additional mechanism besides dual conformal 
symmetry fixes their form at six points and beyond. 

PACS numbers: ll.15.Bt, 11.15.Pg, 11.25.Db, ll.25.Tq, 12.60.Jv 
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I. INTRODUCTION 



Gauge theories play a central role in modern particle physics. How far can we go in 
understanding their properties quantitatively? Long ago 't Hooft suggested that the 
theories simplify dramatically in the so-called planar limit, that of a large number of colors 
N c . In this limit, he suggested that the theory is expressible as a string theory. This idea was 
given a concrete realization by Maldacena's [2] conjecture of the anti-de Sitter/conformal 
field theory (AdS/CFT) duality between weakly-coupled type-IIB string theory on an AdS 
background and the maximally supersymmetric (A/ - = 4) gauge theory at strong coupling. 

The duality requires that the perturbative series for a variety of quantities, including 
scattering amplitudes, sum up to a simple quantity. In order for this to be possible, it would 
appear almost essential for the terms in the perturbative series to be related to each other 
in a simple way. Such a relation is guaranteed for the infrared- singular terms in amplitudes 
by the requirement that physical quantities be infrared-finite. However, the finite terms are 
not required to obey such an iterative relation. The Af = 4 theory is special, as was shown 
by Anastasiou and three of the authors (ABDK) |3|], in that a simple relation does hold 
between the finite terms of the planar one- and two-loop four-gluon amplitudes. (Using the 
Ward identities of Af = 4 sirpersymmetry, this relation extends to four-point amplitudes for 
arbitrary external states 0, By demanding that the amplitude have proper factorization 
as any two external momenta become collinear, this relation was extended to an arbitrary 
number of external legs, valid for the simplest configuration of gluon helicities, the maximally 
helicity violating (MHV) one. Smirnov and two of the authors (BDS) formulated [s] an all- 
orders version of this relation, based on the idea that in the Af = 4 theory the finite terms 
should obey the same (exponential) relations as the divergent terms. This ansatz gave a 
prediction to all orders in the coupling for planar MHV amplitudes. 

In a remarkable paper, Alday and Maldacena [?J suggested a way to compute the 
dimensionally-regulated planar four-point amplitude at strong coupling, using dual string 
theory. Their calculation reproduced the strong-coupling limit of the BDS ansatz. 

The Alday-Maldacena calculation also provided a link between scattering amplitudes at 
strong coupling and a special kind of Wilson loop, one composed of light-like segments in 
a dual coordinate space. Drummond, Korchemsky, and Sokatchev js] found the surprising 
result that at weak coupling the lowest-order (one-loop) contribution to a rectangular Wilson 
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loop with light-like edges is equal to the four-point one-loop amplitude (normalized by the 
tree amplitude, and up to a constant term). Brandhuber, Heslop, and Travaglini [9J showed 
that the equality extends to one between Wilson polygons with n sides and one-loop MHV 
n-point amplitudes for all n. 

Does this equality extend beyond one loop? The perturbative planar amplitudes exhibit 



a "dual conformal invariance" 
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111 ], distinct from the usual conformal invariance of the 



M = 4 theory. Roughly speaking, it corresponds to conformal invariance in momentum 
space. Although the origin of this symmetry is not yet understood, we can use it to guide 
the calculations. Scattering amplitudes with more than four gluons, when normalized by 
the tree amplitude, may contain odd powers of the Levi-Civita tensor contracted with the 
external momenta. These terms flip sign under a parity transformation, which reverses all 
helicities. We refer to them as odd terms, and to the remaining terms as even. One may write 
the four-point amplitude through five loops, and the even part of the five-point amplitude at 
two loops, purely in terms of pseudo- conformal integrals. These are dimensionally- regulated 
integrals whose off-shell continuation is invariant under conformal transformations of dual 
coordinates, whose differences are momenta. 

As shown by Drummond, Henn, Korchemsky, and Sokatchev (DHKS), the dual conformal 
symmetry gives rise 12|, [13( to an anomalous Ward identity in dimensional regularization. 
This Ward identity fixes, to all loop orders, the finite parts of the expectation values of the 
Wilson loops corresponding to four- and five-point scattering amplitudes (up to a constant 
term interpretable as an ultraviolet scale adjustment). The same authors also showed that 
these unique solutions of the anomalous Ward identity coincide with the four- and five- 
point two-loop Wilson loops and amplitudes as computed in perturbation theory (and again 
normalized by the tree amplitude). Given as well the one-loop equivalence between Wilson 
loops and amplitudes 8, 9] and the strong coupling results of Alday and Maldacena 7], we 
may expect that the same equivalence holds for the four-point amplitude and the even part 
of the five-point amplitude, to all orders in perturbation theory; indeed, into the strong- 
coupling regime. The BDS ansatz is in fact a solution to the anomalous Ward identity 12]. 
Because of the uniqueness imposed by the Ward identity, the BDS ansatz too should be 
expected to hold to all orders for four- and five-point amplitudes. It has been checked for 
the four-point amplitude through three loops 0, 0] and for the five-point amplitude through 
two loops fl4. [l5| . 
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What about amplitudes with a larger number of legs? Alday and Maldacena have 
shown 16| that in the limit of a large number of legs, the Wilson loop calculation does 



not agree with the BDS ansatz. This result might imply that the connection between Wil- 
son loops and the amplitudes breaks down; or it might mean that the BDS ansatz breaks 
down for more than five external legs. But at how many legs and how many loops might 
the breakdown occur? 

It is possible to examine the ansatz for consistency in different kinematic limits. Recently, 
the BDS ansatz has been examined in various types of Regge, or high-energy, limits of 



the scattering 



limits 
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he four- and five-point amplitudes appear to be consistent in all such 



191 ]. Indeed, higher-order coefficients in the Regge slope parameter and 



other high-energy quantities can be extracted from such limits. On the other hand, study of 
a particular multi-Regge limit of 2 — > 4 scattering, and also of 3 — > 3 scattering, appears to 



indicate a difficulty with the ansatz for the six-gluon amplitude starting at two loops [19j. 

In order to test the BDS ansatz directly, we have computed the parity-even part of the 
two-loop six-point MHV amplitude in the M = 4 supersymmetric gauge theory. With 



2d 



21], we can 



assistance from the work of Drummond, Henn, Korchemsky, and Sokatchev 
also test the correspondence of the parity-even part with the calculation of a hexagonal 
Wilson loop. Six external legs marks the first appearance of cross ratios invariant under the 
dual conformal transformations; the finite part of the Wilson loop is no longer fixed by the 
anomalous Ward identity, but is determined only up to a function of these cross ratios. Six 
external legs also marks the first appearance of non-MHV amplitudes. The basic Wilson 
loop is insensitive to the helicities of the external gluons; hence it cannot equal these other 
six-point helicity amplitudes, even at one loop. The question of an iterative or exponential 
structure for the non-MHV amplitudes is an interesting one, but we shall not explore it in 
the present paper. 

Another open question, not addressed in this paper, is the behavior of the parity-odd 
part of the two-loop six-gluon MHV amplitude. In the five-gluon case, the parity-odd part 
vanishes in the logarithm of the full amplitude, due to a cancellation between one- and 

rfn 

two- loop terms [14L Ha] . We expect the same cancellation to take place for MHV amplitudes 
with six or more gluons, but this expectation remains to be established. Strictly speaking, 
the parity-odd cancellation is required to establish full correspondence with a Wilson loop, 
which obviously does not change with the reversal of all gluon helicities. 



4 



We perform the calculation using the unitarity-based method, employing a variety of four- 
dimensional and D-dimensional cuts (with D = 4 — 2e) to express the amplitude in terms of 
a selected set of six-point two-loop Feynman integrals. The result may be expressed as a sum 



of pseudo-conforma 
fiveloopsQ, 11, 22, 



loops 
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integrals 10J, in close analogy with the four-point amplitude through 



23 



241 ] and the parity-even part of the five-point amplitude through two 



151 . |25j. There are some additional integrals in the one- and two- loop six-point 
amplitudes, whose pseudo-conformal nature is less clear. Their integrands vanish as D — > 4, 
yet their integrals can be nonvanishing in this limit. However, their one- and two-loop 
contributions conspire to cancel in the logarithm of the amplitude (which is what really is 
needed to test the BDS ansatz) in the limit e — > 0. We then evaluate the integrals using 

n n 

the packages AMBRE [26] and MB [27[ and compute the amplitude numerically at a variety of 
kinematic points. The structure of the infrared singularities is known |28L l29| , and agrees 
with the pole terms in our expression. The finite remainders are tested numerically against 
the BDS ansatz, and against values for the corresponding Wilson loop 21]. 



The paper is organized as follows. In section [TTJ we review the ABDK and BDS ansatze, 
the structure of scattering amplitudes at strong coupling, the difficulty that appears as the 
number of external legs becomes large, and dual conformal invariance. In section IITT1 we give 
the integrand, and outline its calculation via the unitarity method. In section HVl we present 
our results for the amplitude and compare these to the results of the hexagonal Wilson loop 
calculation. Section |V] gives properties that the remainder function must satisfy. We give 
our conclusions and summarize open problems in section [VTl The appendices contain results 
for the integrals. 

II. REVIEW 

A. ABDK/BDS Ansatz 

In an SU(N C ) gauge theory the leading-color contributions to the L-loop gauge-theory 
n-point amplitudes can be written as 



A™ = ^- 2 a L ^Tr(T^)T^)...T a ^))4 L )(p(l),p(2),...,pH) 



(2.1) 
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where 

A 

Here A = g 2 N c is the 't Hooft parameter, g is the Yang-Mills coupling, and 7 is Euler's 



a = (4ne^y— . (2.2) 



constant. The sum is over non-cyclic permutations of the external legs. We have suppressed 
the momenta and helicities fc, and Aj, leaving only the index i as a label. This decomposition 
holds for any amplitude when all particles are in the adjoint representation. We will find it 
convenient to scale out the tree amplitude, defining 

M^\e)=A^/Af. (2.3) 

In planar M = 4 supersymmetric gauge theory, amplitudes computed to date satisfy 
an iteration relation. At two loops, the iteration conjecture expresses n-point amplitudes 
entirely in terms of one-loop amplitudes and a set of constants j^]. For two-loop MHV 
amplitudes the ABDK conjecture reads 

Mf(e) = ~(M«(e)) 2 + f®(e) M«(2e) + C* a > + {e) , (2.4) 

where 

f i2 \e) = -(Ca + C 3 e + <U<? + ■■•), ^ = . (2.5) 

The form (12. 4p was based on explicit computations of both the four-point two-loop ampli- 

3j, which control the behavior of the amplitudes 



tude [23j, [3Cj and of the splitting amplitudes 
as two momenta become collinear. The collinear splitting amplitude has an iterative prop- 
erty, analogous to eq. 02.4)) . which guarantees that the ansatz 02.4)) for MHV amplitudes 
has the correct collinear limits for any n. In addition, MHV amplitudes have no poles as 
a single multi-particle kinematic invariant vanishes, (hi + • • • + kj) 2 — > for j > (i + 1) 
mod n; so factorization in such channels is satisfied trivially by eg. 02 -41). For the five-point 



amplitude, eq. 02 .4p has also been confirmed by direct calculation [14j, [l5j. Any violation of 
eq. 02.4)) beyond five external legs must necessarily be expressed by a function which van- 
ishes as pairs of color-adjacent momenta become collinear. We shall see in section M that 
such a remainder function is, however, detectable in the triple-collinear limit in which three 
color-adjacent momenta become collinear. This limit can first be achieved for n = 6. 

The iterative structure in eq. 02.4)) . together with the exponential nature of infrared 



divergences 



28 



29], suggest that an all-orders resummation should be possible. In ref. 



16] 
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the three-loop generalization for n = 4 was found by direct calculation, guiding the all-loop 
order BDS proposal, 



In M n = a ' (f (l) ( £ ) M n 1} fa) + C (l) +0{t 



1=1 



where 



(2.6) 



(2.7) 



L=0 



is the resummed all-loop amplitude. The quantity Mn\lt) is the dimensionally-regulated 
one-loop amplitude, with the tree scaled out according to eq. ( 12. 3D . and with e — >■ /e. Each 
f^'(e) is a three-term series in e, beginning at C(e°), 



(2.8) 



The constant /q is the planar cusp anomalous dimension 31], = 7 



In order to test the ABDK/BDS ansatz, it is convenient to define an /-loop remainder 
function B$ to be the difference between the actual /-loop rescaled amplitude M„ and the 
ABDK/BDS prediction for it, in the limit e — > 0. This function is finite as e — > 0, because 
the BDS ansatz has all the correct infrared singularities. It is only defined in the limit e — > 
because the two-loop ansatz does not hold beyond O(e ), even for n = 4. For example, at 
two loops the remainder function is defined by, 



R {2) = lim 



M%\e) - (\{M£\e)f + f®(e) M«(2e) + C<») 



(2.9) 



Notice that the combination Mn^ — \/2{Mn) 2 appearing in eq. (12. 9p is the order a 2 term 
in the logarithm of the amplitude ( 12.61) . 

oop MHV amplitudes entering the ABDK/BDS ansatz were computed some 



The one- 



time ago |32|], with the result 



+ Fi 1 )(O) + 0(6) 



(2.10) 



11] . The "0" argument in the finite part 



where we use the normalizations of refs 
Fn\fy signifies that we have taken e — > 0. These terms have the form, 



n 



'2.11 



i=l 



where 



9n,i 



Ln/2J-1 



r=2 



I*^=£ ) ln( ~ g(m) - ( ^M + Am + L n>i + ^2 , (2.12) 



in which |_^J is the greatest integer less than or equal to x. Here Si...j = (&$ + ••• + fcj) 2 are 
the momentum invariants. (All indices are understood to be modn.) The form of D n ^ and 
L n i depends upon whether n is odd or even. For the even case (n = 2m) these quantities 
are given by 

m— 2 ✓ 
£>2m,i = ~ ^ Li 2 ( 1 



r=2 



S»...(i+r— (i+r) \ 1 T . / s i---(i+m-2) s (i-l)---(i+m-l) 



s i— (i+r) s (i— 1)— (i+r— 1) 



Li 2 1 1 

- V Sj...(j +m _i)S(j_i)...(j +m _2) 



r 1 i„2 / s i--(i+m-l) 

^2m,i — T m I 

4 \-S(j + i)...(j +m ) 

In the odd case (n = 2m + 1), we have, 



(2.13) 



m— 1 



2m+l,i 



EM 1 



r=2 



li—(i+r-l) s (i-l) — (i+r) 



2 \ "Si---(j+m) 



s i---(i+r) s (i-l)---(i+r-l) . 

-S(i+l)-(i+m) 



"■S(j-l)...(j+m-l) . 

For n = 4 the above formula does not hold; in that case the finite part is simply 

^ (1) (0) = ^ln 2 (3+4C2. 



(2.14) 



(2.15) 



To make contact with the string theory literature [3, H, |33j we define 



/(A) =4 a'/, 



j/(0 
o J 



</(A) = 2£^/< 



a .(I) 



i 00 i 



Z=l 



Z=2 



(=2 



In terms of these functions, the BDS ansatz (12. 6p may be written as, 

\nM n = Div n +^F«(0) + nk(\) + C(X) . 



(2.16) 



(2.17) 



The infrared-divergent part is 

n r i 

^ = -£[^/ ( - 2) (fZ 

where 



j=i 



A/z 2e 



1 ' ^ 2e 



^)>(A) = /(A) 



4e y 



A^)^ 1) (A)=^(A), 



A/i 



~ s i,i+lj 



(2.18) 



(2.19) 



and n 
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m 



Ane "<fi 2 . The first few orders of both the weak 



29 



34 



351 ] and strong 



36, 



381 ] coupling expansion for /(A) have been computed, with the result 



/(A) = 


A 




11 A 2 


/ 73 


2tt 2 


11520 


V 1290240 


/(A) = 


Vx 

IT 


f 31n2 
V 1 v/A 


K 

-x + " 


•)■ 



+ 



CI 



5127T 6 



A^ + 



A -> 



oo 



(2.20) 
(2.21) 



where K = 



(-i)" 



, „ „ , l7 r _= 0.9159656. . . is the Catalan constant 1 . The function q(X) has 

' n >° (2n+l) 2 i-i I — I n 

also been computed through four loops [g, |40j, but is less well known at strong coupling |7|] 2 , 



fl'(A) = -C 3 (^) 2 + ^(6C5 + 5C 2 C3] 
g(X) = (l-ln2)^ + ---, 



A 



(77.56 ±0.02) 



A 



A 
A - 



•0, (2.22) 
oo. (2.23) 



Beisert, Eden and Staudacher (BES) proposed [33|, |4l| a striking integral equation giving 
the cusp anomalous dimension /(A) for all values of the coupling. This integral equation has 



35| and strong coupling [37, 



39 



42|, 



passed a number of stringent tests at both weak ll|, 
and is therefore very likely the correct expression. 

With the cusp anomalous dimension known, the BDS ansatz (12.61) predicts the MHV 
amplitudes for all values of the coupling, up to the undetermined functions g(A), k(X) and 
C(A), which are independent of the kinematics. The ansatz has been checked through three 

n n 

loops at four points J3L [6J . Integral representations of the four-point amplitude have also 
been given at four [llj and five loops 24|, though these expressions have not yet been 
integrated at O(e ) to yield explicit functions of the external momenta. If one assumes the 
dual conformal invariance mentioned in the Introduction, which we shall discuss at greater 
length below, then the form of the finite parts in the four- and five-point amplitudes are 
fixed, which provides another way of arriving at the ansatz (12. 6p for n = 4 or 5. However, 
beyond five external legs, the assumption of dual conformal invariance does not suffice to 
fix the functional form of the finite parts of MHV amplitudes. 



1 The third term in eg <j^.2ip was first found in the expansion of an integral equation 33( in ref. 38 1 . 

2 At the next order [38( in the strong-coupling expansion of g(X), one encounters difficulties with the 
closed-string version of dimensional regularization used in ref. 
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B. Scattering Amplitudes at Strong Coupling 



Alday and Maldacena have proposed a very interesting way to compute color-ordered 
planar scattering amplitudes at strong coupling, using the AdS/CFT correspondence 
They argue that the leading dependence of any amplitude on the coupling has the form 
exp(— ^A), where A is the regularized area of a special surface in AdS, whose definition 
and properties are reviewed below. Their result reproduces the BDS ansatz for the four-point 
amplitude at strong coupling. 



Long ago, Gross and Mende |43j showed that in string theory in flat space-time, scattering 
amplitudes can be computed in the high-energy and fixed-angle regime using a semiclassical 
approach. Alday and Maldacena noticed that, thanks to the properties of anti-de Sitter 
space, a semiclassical calculation suffices to compute the strong-coupling limit of amplitudes 
in the M = 4 theory, for any energy. Glossing over technical details, scattering amplitudes 
are given by a saddle-point approximation to the world-sheet partition function with certain 
vertex operator insertions. 

Two-dimensional duality transformations map a vertex operator of momentum to a 
null (light-like) segment pointing along the direction of the momentum of the corresponding 
gluon in a dual AdS space. The endpoint coordinates y of each null segment obey 

Ayf = 2irk? . (2.24) 

Momentum conservation then implies that these segments form a closed polygon in these 
dual coordinates. 

Alday and Maldacena do not compute the prefactor of exp(— ^-A), which is subleading 
in the strong-coupling expansion but must contain all the dependence on the polarizations 
of the external particles. MHV scattering amplitudes in the Af = 4 theory are special 
because supersymmetry Ward identities 4j, |5( imply that they are all identical up to simple 
spinor product factors. Equivalently, MHV amplitudes can be written as a product of the 
tree amplitude, and an additional factor dependent only on the momentum invariants, and 
not on the polarizations. In the MHV case, it is natural to identify the tree amplitude 



with the prefactor. (A similar proposal has also been made for the non-MHV case 44j, 



although it is not clear how it can be consistent with the intricate structure of the one-loop 



amplitudes 45|, |46(.) Then the additional factor is given, to leading-order in the strong- 



coupling expansion, by the saddle-point approximation of a suitably restricted world-sheet 
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partition function. The world sheets must be restricted to unpunctured surfaces whose 
boundary is the closed polygon of null dual segments. That is, up to some technical details, 
the factor is just the expectation value of a null Wilson loop in the dual space. 

Alday and Maldacena computed the leading term in the strong-coupling expansion of the 
four-gluon amplitude by building an explicit solution out of cusp solutions written down 



earlier by Kruczenski 
equations 



361 ] . The solution can also be found by solving the sigma- model 



47j once the Virasoro constraints are imposed [48j. The minimal- surface ap- 



proach has not yet yielded complete expressions for higher-point amplitudes, because of 
difficulties in solving the minimal surface conditions with the proper boundary conditions, 
although there has been some progress 47j, |49|. However, Buchbinder 50| has shown that 
the infrared- divergent terms can be obtained using this approach, and that (as expected) 
they are consistent with the known exponentiation in gauge theories. Komargodski has 



argued 



5l| that the collinear-splitting amplitude obtained from this approach is consistent 



with the known perturbative splitting amplitudes. 

While the original arguments have been formulated at strong coupling, the work of Al- 
day and Maldacena also inspired the discovery of a surprising connection between MHV 
amplitudes and Wilson loops at weak coupling at one loop (si, jj], for the simplest two- loop 
amplitudes Q , and perhaps also at higher orders. 



C. Trouble at Large n 

In the limit of a large number of external legs, Alday and Maldacena [3] argued that 
the situation simplifies because one may approximate the null Wilson loop by a smooth one. 
They considered a rectangular Wilson loop with space-like edges of length L and width T as 
the approximation of a null Wilson loop zig-zagging around the rectangular one. Each edge 
corresponds to many gluons moving in one direction, alternating with many gluons moving 
in an opposite direction. The expectation value of the Wilson loop is clearly divergent as 
the area is scaled to infinity. However, one may focus on the scale-invariant part, which is 
proportional to T/L. Moreover, one may consider the kinematic configuration corresponding 
to T/L ^> 1. In this case one may ignore the contribution from the sides of length L to 
leading order. Thus one may further approximate the Wilson loop by two parallel lines of 
length T at distance L from each other. In other words, the dominant T/L-dependent part 



11 



of the world-sheet area is essentially T times the heavy-quark potential, computed at strong 
coupling [521 ] : 

HW) = ^ ^4 T ! A » 1 . (2.25) 
r (l) 

Since this expression is finite as T and L are scaled to infinity, one should compare it with 
the T/L-dependent terms in the logarithm of the finite part of the BDS ansatz. 

Alday and Maldacena worked out the behavior of the BDS ansatz in this limit 
making use of the all-n one-loop relation between MHV amplitudes and Wilson loops 
and the known value of the rectangular Wilson loop at one loop. They obtained, at strong 
coupling and for T 3> L, 

4 4 L 

This result differs from eq. (12.251) . indicating that the BDS ansatz is incomplete in this limit. 
It is interesting to note that the two formulae differ by less than 10%, hinting that it may 
be possible to systematically correct the BDS ansatz. 

While this difference appeared in a particular (and somewhat singular) kinematic con- 
figuration with a very large number of legs, the fact that it depends on T/L makes clear 
that it arises from a nontrivial function of momenta. As we discuss below, dual conformal 
symmetry suggests that the first place such a function can occur is at six points. 




D. Pseudo-Conformal Integrals and Dual Conformal Invariance 



We now turn to a review of the observed restrictions that dual conformal invariance places 



on the integrals appearing in the amplitudes [8|, [10|, llll, |24| . This mysterious symmetry plays 
an important role in our story. M = 4 super- Yang-Mills theory is a conformal field theory 
at the quantum level; conformal invariance may be observed in correlation functions of 
operators of definite (anomalous) dimension. However, the constraints it imposes on on- 
shell scattering amplitudes are obscured, both by the need for an infrared regulator and by 
anomalies analogous to the holomorphic anomaly of collinear operators 53]. 

By inspec ting the known results for the one-, two- and three- loop four-gluon ampli- 



tudes 



231 ]. Drummond, Henn, Sokatchev and Smirnov 10J observed that after con 



tinuing the external momenta of the integrals off shell (in a sense we will describe below) 



12 



and then taking D = 4, they exhibit an 5*0(2,4) dual conformal symmetry. This symmetry 
is distinct from the four- dimensional position-space conformal grou p. Interestingly it holds 
individually for each contributing integral through five loops 0, 11, 24]. A discussion of the 



consequences of dual conformal symmetry for complete amplitudes, instead of individual 



integrals, may be found in refs. 12J, Il3 ] 



The origin of dual conformal symmetry remains obscure, and its broad validity for am- 
plitudes remains to be proven. 3 Indeed, it is not clear that dual conformal symmetry holds 
for all contributions to scattering amplitudes. For example, at the five-point level, the dual 
conformal properties of the parity-odd pieces are not apparent 151 ] . (It is conceivable that 
they might be re-expressed as a linear combination of integrals with manifest properties, 
although it seems unlikely, because the single available parity-odd contraction of the Levi- 
Civita tensor at the five-point level does not transform homogeneously under the inversion 
of dual variables described below.) At least in this case, the parity-odd pieces do not enter 
the remainder function R§ defined in eq. (12.91) . This feature suggests that, as long as we 
subtract appropriate lower- loop contributions (to be specified further below), we may be able 

(2) 

to identify simple conformal properties for the finite remainders R n of MHV amplitudes to 
all loop orders. In this vein, an anomalous conformal Ward identity has been proven for the 



finite remainders of Wilson loops jl3j ]. and later found at strong coupling [511 ] . The BDS 
ansatz was shown to obey the same Ward identity, which was then proposed to hold for 
MHV amplitudes [lj. 



Dual conformal symmetry is most transparent in terms of dual variables Xi which are 
related to the gluon momenta in a way analogous to eq. (12.241) . 

ki = a^j+i — Xi , (2.27) 

and similarly for the loop momenta. Formally, in an integral we may identify the variables 
Xi as the positions of the vertices of a dual graph. In this construction the momentum con- 
servation constraint is replaced by an invariance under uniform shifts of the dual coordinates 
Xi. Since the dual variables are unconstrained and the parametrization (12.271) automatically 



3 Dual conformal symmetry was introduced in the context of multi-loop ladder integrals [54|; it has also 
cropped up in the two-dimensional theory of Reggeon interactions, again in the planar limit \§l 
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x 3 



X 2 tH 



•r,5 



(a) 




k 



3 *4 k 



X 2 



If X7 




,x 5 



(d) 



FIG. 1: Examples of pseudo-conformal integrals. Points Xi label the vertices of the dual graph, 
a solid line connecting two points X{ and Xj corresponds to a factor of 1/xfj, while a dashed line 
corresponds to a factor of xfj. An integral is pseudo-conformal if the difference between the number 
of solid lines and dashed lines at a vertex equals 4 at the internal vertices and zero at the external 
vertices. Graphs (b), (c) and (d) show that the integrals appearing in the even part of the five-point 
two-loop amplitude are pseudo-conformal. (In fact, (d) only appears in the odd part.) 

satisfies momentum conservation, one may define an inversion operator 



xj 



(2.28) 



One may also define conformal boost transformations of the dual variables x; they are 
generated by 



i ' 



^2.29) 



Invariance under I implies invariance under K, because the conformal boosts are generated 
by two inversions, with an infinitesimal translation of x in between. 

The principal conformal-invariance constraints on integrals constructed from the invari- 
ants xfj are exposed by performing the inversion I. Since dimensional regularization breaks 
the dual conformal invariance, for the purposes of exposing the symmetry we adopt a differ- 
ent infrared regularization of the integrals. We take the external legs off shell, letting kf ^ 0, 
instead of using a dimensional regulator as in the rest of the paper. Under the inversion, 
the Mandelstam invariants and the integration measure transform as 



X; 



X 



'.I 



1 2 ' 

xfxj 



d X5 



d A x^ 



x 



2\4 



d 4 XQ 



d A x§ 



x, 



2\4 



^2.30) 



Simple dua 



refs 



. [lO 
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diagrams may be used to identify all the integrals that are invariant under /; see 



24j for a detailed discussion. We call such integrals "pseudo-conformal" because 
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jlm j2me jhex 

FIG. 2: The three independent integrals which contribute to the six-particle amplitude at one loop. 
The external momenta are labeled clockwise, beginning with k\ which is denoted by an arrow. The 
Hp in the hexagon integral indicates that a factor of the square of the (— 2e)-dimensional components 
of the loop momentum is to be inserted in the numerator of the integrand. 

in the amplitudes we actually use the dimensionally-regulated versions of the integrals, which 
breaks the conformal invariance but guarantees the gauge-invariance of the amplitudes. 

While pseudo-conformal integrals were initially identified in four-point scattering ampli- 
tudes, they apparently play a more general role in the planar J\f = 4 theory; they also appear 



in n-point one-loop MHV amplitudes 
amplitud e fl4 
one loop 



321 ] , as well as the even part of the two- loop five-point 
Il5j ] (see fig. [I]). Their relevance for non-MHV amplitudes, for example at 
45] , remains to be clarified. In this paper we show that for the two- loop six-point 
MHV amplitude, they again appear in the even part though some additional integrals also 
appear. (We shall not compute the odd part in this paper.) 

More generally, conformal invariance implies that, as for two-dimensional conformal field 
theories, conformally invariant quantities can depend only on the cross ratios, 

X ■ cc^ 

= -jf-2- ' ( 2 - 31 ) 

x ik x jl 

where i,j, k and I are any four vertices of the dual graph. In scattering amplitudes we 
must exclude cross ratios that vanish or diverge on account of an x^ vanishing due to the 
on-shell conditions, xj i+l = kf — > 0. Consequently, the only cross ratios that may appear 
in a conformally invariant part of an on-shell amplitude are those with \i — j\ > 2 modn 
for all i,j — 1, . . . , n; otherwise would correspond to a massless momentum. We remark 
that there are no conformal cross ratios for n = 4, 5, as it is not possible to form a cross 
ratio without encountering a vanishing xf„- [1 
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III. SIX-POINT INTEGRAND 



Before turning to our calculation of the two-loop amplitude, we present the one-loop six- 
point amplitude (5j. We drop the parity-odd terms, which are proportional to the Levi-Civita 
tensor and are also 0(e). We split the even part into two pieces, 

M (l),D=4-2e (e) = M (D^=4 (e) + M (D^ (e) j (3 1} 



where 

„s r, , i _ r . 1 1 

(3.2) 



M^ D =\e) = -\ £ 



'6 4 

b perms. 

and 



S45S 56 / lm (e) + i(si 23 S 3 45 - Si 2 S 45 )/ 2me ( 



M 6 (1) ' M (e) = -- tr[123456]/ hox (e) . (3.3) 

The first piece contains scalar box integrals, which are constructible solely from cuts in 
D = 4. The hexagon integral in the second piece contains a numerator factor of fi 2 , as 
indicated in fig. [2, which can only be detected by cuts in which the cut loop momentum p 
has a nonvanishing (— 2e)-dimensional component p 2 = p? 4 , — fi 2 , so that /J, 2 is positive. 
(The overall scale fi arising from dimensional regularization, e.g. in eq. (12.101) . should not 
be confused with /i p .) 

The three integrals appearing in the one-loop amplitude are defined in fig. [2, using the 
convention that each loop momentum integral is normalized according to 



m 



- D ' 2 e^ d D p. (3.4) 



The coefficient of the hexagon integral involves the quantity 

tr[123456] = tr^^hhhh] = S123S234S345 - s 6 is 34 si23 - S12S45S234 - S23S56S345 • (3.5) 

The sum in eq. (13 .2p runs over the six cyclic permutations of the external momenta k{. There 
is a symmetry factor of 1/2 in front of I 2me to correct for a double count in this case - 
there are only three separate J 2me integrals. 



A. Construction of the Integrand from Unitarity Cuts 



We computed the two- loop six-point MHV amplitude using the unitarity method 32, 



45 



561 ]. We know on general grounds that the amplitude can be written in terms of Feynman 



integrals multiplied by rational coefficients. 



16 




y 



(6) 




\ 








\ 


/ \ 






























/ 


\ 




\ 









(2) 



(7) 
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x(p + h) 2 {q + h) 2 x(p + k 3 ) 2 (q + k 2 ) 2 

(12) (13) (14) (15) 

FIG. 3: The 15 independent integrals which contribute to the even part of the six-particle amplitude 
at two loops. The external momenta are labeled clockwise with k\ denoted by an arrow. Integrals 
(8)-(15) are defined to include the indicated numerator factors involving the loop momenta. In 
the last two integrals, fi p denotes the (— 2e)-dimensional component of the loop momentum p. 



In principle, we should consider all two-loop six-point integrals. We can reduce this 
rather large set of integrals using the no-triangle constraint, which states that one can find a 
representation for the amplitude in which no integral with a triangle (or bubble) subintegral 
appears. Using generalized unitarity, we may establish this constraint for two classes of 
potential triangle contributions: those which may be excised from a multi-loop topology by 
cutting only gluon lines, and (using sup ersymmetry Ward identities) those which may be 
excised by cutting at most five legs llj. We shall assume that the no-triangle constraint 
holds as well for the other contributions. 

We expect, based on previous calculations of four-point amplitudes at two, three, and 
four loops (as well as the consistency of the five-loop construction) and of the five-point 
amplitude at two loops, that essentially only pseudo-conformal integrals will be required for 
the six-point amplitude. The basic topologies for the integrals that appear in the planar six- 
point MHV amplitude are shown in figs. [3] and HI Of these, all are pseudo-conformal, except 
for (14) and (15), whose integrands vanish as D — > 4, yet integral (15) is nonvanishing in 
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V 




(16) 



(17) 



x(p + ki) 2 
(18) 

FIG. 4: The three independent two-loop diagrams which can be made pseudo-conformal by includ- 
ing appropriate numerators but which do not contribute to the amplitude (i.e., they enter with 
zero coefficient). 



this limit. The complete set of inequivalent pseudo-conformal integrals (including numerator 
factors) is presented in fig. [7J The no-triangle constraint would also follow from the stronger 
assumption of pseudo-conformality, were we to make it. We will not; our calculation can be 
seen as a test of it instead. The appearance of integrals (14) and (15) shows that additional 
pieces, unobvious from naive considerations of dual conformal invariance, enter. 

Given the no-triangle constraint, iterated two-particle cuts suffice to determine the full 
integrand. While the four- dimensional cuts do not suffice to determine the integrand com- 
pletely, they can determine most terms. The four-dimensional double two-particle cuts, 
shown in fig. \5\ are particularly convenient to compute because they can be built out of 
MHV tree amplitudes. 

As example, consider the helicity assignment (1~, 2~, 3 + , 4 + , 5 + , 6 + ), and compute cut 
a. Because of a Ward identity for J\f = 4 supersymmetry 0, once we have divided by 
the tree amplitude, the expression is in fact independent of the placement of the negative 
helicities. The labeling of the external legs is shown in fig. |6j 

The product of the three tree amplitudes corresponding to cut a is 



(12Y 



x i- 



(H 2 ) nor 



X I- 



(hhY 



(2 3) (3 h) (h l 2 ) (l 2 1) ((-JO HO) (HO Hs)> (Hs) H)) (h 4) (4 5) (5 6) (6 1 3 ) ' 

(3.6) 

After spinor simplifications, dividing by the tree amplitude, and rationalizing denominators 
to Lorentz products, we find that the numerator can be written as follows, 



tr+ [I2 $i $6/3/4 ^4 $3/1] tr + [/2/3/4/1] 



(3.7) 



where tr + [- ■ ■] = |tr[(l +75) - • •]■ 

Upon expanding the traces, we find both even and odd terms. The odd terms contain a 
factor e(a, b, c, d) = e^ upa a^b u c p d a , whose origin lies in the presence of the 75 matrix inside 
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1 1 II 1 1 

(a) (b) 

I I I I I I 



(c) 



I I I 



FIG. 5: The double two-particle cuts used to determine the integrand. 

h ■ A II < < -4^ ;'— A ' A'4 

^2 < — H H H — > fa 

A'l — >- ; — — /- A'g 

FIG. 6: Labeling of momenta for cut a. 

the traces. (The product of two epsilon tensors would yield an even term, but only the 
longer trace here can actually produce an epsilon tensor, as e(7i, l 2 , h, h) vanishes because 
of momentum conservation.) We ignore the odd terms in our calculation. 

In order to identify the coefficients of the integrals in fig. [31 we should use momentum 
conservation to re-express all Lorentz invariants in terms of independent invariants. The 
required simplifications can be done analytically, but in some cases (for example cut d) it is 
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easier to do them numerically, by matching to a target expression. 

Doing so, we obtain for the final result of cut a, in the (3, 4, 5, 6, 1, 2) permutation with 
respect to fig. [3l 

S 123 S 34S61 ~ S 123 S 234'S345 + 5i23 s 2345l2545 + S123S345S23S56 



2 



2 2 

_|_ s 123 s 345 — <Sl23"Sl2<S45 _|_ S 123 S 234 — 5i23523 s 56 



(k 3 + h) 2 (h + h) 2 (h + h)2 (h + i 2 ) 2 (i2 + h) 2 (h + k) 2 

I g 123 g 34 S 123 S 61 

/ T T\0/I . 7 \ O / 7 J \ O ~ ' 



+ 



(A* + ^i) 2 (^2 + h) 2 (h + k) 2 (h + h) 2 (l2 + h) 2 (h + h) 2 

5123^12523(^6 ~ hf Sl23Sl2«23 (h - h) 2 



(h + h) 2 (h + h) 2 (h + h) 2 (h + h) 2 (h + l 2 ) 2 (h + h) 2 (h + h) 2 (h + U 



\2 



-§123545556(^3 ~ U) 2 5i 2 35 4 5556 (fa ~ h 

/■» 1 \ n / t 1 \ n / 1 1 \ o / 1 7 \ n ^ 



(A* + /i) 2 (^2 + h) 2 (h + h) 2 (h + h) 2 (fa + / 2 ) 2 (/ 2 + h) 2 (fa + / 4 ) 2 (fc 6 + h) 2 

1 

+ (fa + h) 2 {h + /i) 2 (/ 2 + h) 2 (fa + h) 2 (fa + h) 2 



'123 



5 6 i(^ 3 - h) 2 (fa - k) 2 - s 2 l23 s u {fa - / 3 ) 2 (^ 6 - kf 



+ 5123(51235234 - S23S56X&3 - h) 2 (fa ~ h) 2 

+ 5i23(5l235345 — 5i2S45)(fci — Z 3 ) 2 (fc 4 — li) 2 



(3.8) 



From this expression we can read off the coefficients of every integral in figs. [3] and H] that 
is nonvanishing on cut a for four- dimensional values of the cut momenta. These integrals 
are, in order of their appearance in eq. (EED, /W, (twice), J( 4 > (twice), /( 9 ) (four times), 
f( 12 ) and i^ 13 ** (each twice). For example, from the second term, up to a symmetry factor, 
we can simply read off the coefficient of 1^ to be 5123(51235345 — S12S45), or equivalently 
5234(51235234— 523S 56 ) in the labeling of fig. [3) The numerator in the second term on the fourth 
line produces the required loop-momentum dependent factor for a reflected version Of I®, 
corresponding to the (4, 3, 2, 1, 6, 5) permutation. The coefficient is 5i 2 35i 2 5 2 3, or S234534S23 
in the figure's labeling. The coefficients of the remaining integrals can be determined by the 
other cuts in fig. [51 

At one loop in any supersymmetric theory, the improved ultraviolet power-counting en- 
sures that any rational terms in an amplitude are linked to terms with branch cuts. That is, 
all terms in an amplitude can be determined solely from the standard integral basis, which 
can be detected in four-dimensional cuts. Beyond one loop, four-dimensional cuts no longer 
suffice for J\f = 1 supersymmetric theories 53]. In the J\f = 4 theory, four-point amplitudes 
through five loops are determined solely by their four-dimensional cuts. The same is true for 
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the even terms in the five-point amplitude. (It is no longer true for the odd terms, but in any 
case we are ignoring the corresponding terms in the six-point amplitude.) However, there 
is no proof other than explicit computation of this observation. Accordingly, we cannot be 
certain that four-dimensional cuts will suffice for our calculation. 

Indeed, the hexagon integral in fig. [2] will contribute terms of 0(e) to the one-loop am- 
plitude. In the term in the iteration relation (12.41) in which the one-loop amplitude appears 
squared, the product of such terms with singular terms in M^\e) survives to give 0(e~ 1 ) and 
finite contributions. We will see in section IIV Al that such contributions are offset by those 
coming from the last integral (15) from fig. O induced by the (— 2e)-dimensional components 
of the loop momentum. These contributions must be computed using D-dimensional cuts, 
either making use of prior computations Q], or by direct computation. We have computed 
two cuts, corresponding to fig.[5](a) and (c), using D-dimensional cuts. These cuts determine 
the coefficients of integrals (14) and (15), respectively, in fig. [3j The calculations were done 



using the same approach used in refs. 3|, Il5j . While we can no longer use standard helicity 
states for the computation, we can take advantage of the equivalence between the TV = 4 
theory and ten-dimensional M = 1 super- Yang-Mills theory compactified on a torus. We 
compute the cuts with the spin algebra performed in the ten-dimensional theory, keeping 
loop momenta in D dimensions. (External momenta can be taken to be four-dimensional.) 
The ten-dimensional gluon corresponds to a four-dimensional gluon and six real scalar de- 
grees of freedom, while the ten-dimensional Major ana- Weyl fermions correspond to four 
flavors of gluinos. 



B. Presentation of the Integrand 

By analyzing the cuts outlined in the previous section, we find the complete expression 
for the parity-even part of the two-loop six-particle amplitude to be 



M ( 2 ),D=4-2e (e) 



(2),D=4 



(i 



e) + M, 



(2),M/ 



(3.9) 



where 



- E 

16 ^ 



12 perms. 



+ c 6 / (6) ( e ) + \c7l {7 \e) + \c,I { *\e) + c 9 I®(e) + c 10 I^(e) 
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+ c n I {U \e) + ±c 12 I^\e)+ 1 -c 13 I^\e) 



(3.10) 



and 




1 



(3.11) 



16 



12 perms. L 



involving the 15 independent integrals shown in fig. [3j As with the one-loop ampli- 
tude (13 .11) . we have separated the integrals (1)— (13), which are constructible solely from 
four- dimensional cuts, from integrals (14) and (15), in which the (— 2e)-dimensional com- 
ponents of loop momenta fi p and fi q are explicitly present. (The D = 4 superscript in 
Mg 2 ' ) ' D_4 (e) refers to the lack of \i v and fi q terms in the numerators of the integrands; the 
argument e indicates the dependence of the integrals on the dimensional regularization pa- 
rameter.) Of the 1/16 overall normalization in eqs. (13.101) and (13. lip . 1/4 is due to our 
choice of normalization in eq. (12. II) and the other factor of 1/4 emerges from the calculation 
of the unitarity cuts. 

The sum runs over the 12 cyclic and reflection permutations of external legs, 



(1, 2, 3, 4, 5, 6), (2, 3, 4, 5, 6, 1), (3, 4, 5, 6, 1, 2), (4, 5, 6, 1, 2, 3), 
(5, 6, 1, 2, 3, 4), (6, 1, 2, 3, 4, 5), (6, 5, 4, 3, 2, 1), (1, 6, 5, 4, 3, 2), 
(2,1,6,5,4,3), (3,2,1,6,5,4), (4,3,2,1,6,5), (5,4,3,2,1,6). (3.12) 



The numerical coefficients in each term of eqs. (I3.10p and (13.111) are symmetry factors to 
remove double counts in the permutation sum. The coefficients for the (1,2,3,4,5,6) per- 
mutation are, 



c 3 — $234:( s 123 s 234 — S23 s 56) ' 
2 

C4 — Si2S 234 , 

c 5 = ^34(51235234 — 2s 2 3 s 56) j 

Cq = — S12S23S234, 

Cj = 2S123S234S345 — 4S61S34S123 — S12S45S234 — ■S23S56S345 , 

c 8 = 2s6l(s234-5345 — S61S34) , 

c 9 — S23S34S234 , 



s 61 s 34'5l23' s 345 + ■S;L2' S 45' S 234'S345 S 345 ( s 23 s 56 — S123S234) 
2S12S23 > 
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ClO — s 23 (2s 6 i 

Cll = S12S23S234, 

Cl2 = S 345(52345345 — S61S34) , 

c 13 = — S 345 S 56 ) 

Cl4 = — 2^345 (S123S234S345 — S61S34S123 — S12S45S234 — S23S56S345) j 

C15 = 2s 61 (si23S234S345 — S61S34S123 — S 12-S45 s 234 — S23S56S345) . (3.13) 

Equation (13.91) omits all odd terms proportional to the Levi-Civita tensor. Because neither 
the strong-coupling string theory calculations nor the Wilson loop calculations are sensitive 
to such terms, we may drop them without affecting comparisons to either. 

Finally we should note that while the set of integrals in figs. |3] and H] is convenient 
for calculation, non-trivial relations between them and other integrals may allow for other 
equivalent representations of the amplitude (13.91) in which the form of the coefficients (13 . 13[) 
is substantially altered. 



C. Dual Conformal Structure of the Integrand 

As mentioned above, we expect that planar amplitudes should manifest dual conformal 
symmetry. An even stronger statement, that an amplitude can be expressed as a linear 
combination of integrals, each exhibiting manifest pseudo-conformal invariance, has been 
observed to hold for the four-particle amplitude through five loops, and for the even part of 

the five-particle amplitude at two loops. 

(2) 

As with the odd part of M5 , the dual conformal properties of the two integrals in 
Mg 2 ^, i^ 14 ) and are not apparent. It is possible that they can be re-expressed as 

a linear combination of integrals with manifest properties. However, we will argue in the 
next section that they cancel against the analogous one-loop piece Mg , in the remainder 
function ( 12.91) for the ABDK/BDS ansatz. Hence we focus on the surviving piece M^' D ~ A . 
We find that this piece can indeed be written as a linear combination of the 26 independent 
pseudo-conformal integrals, as exhibited in fig. [7J 



The integrals appearing in the four-point amplitude have a number of striking 



partly explained by heuristics such as the rung rule [23] and box substitution rule [24j . Three 



'eatures, 



interesting features are observed in the four-point amplitude through five loops: 



All pseudo-conformal integrals appear with relative weights of ±1 or |3j, QjJ, |23|, \2A\ 

23 




FIG. 7: The 26 different integrals which are allowed, by the hypothesis of dual conformal symmetry, 
to contribute to the amplitude Mg 2 "*' D_4 . Beneath each diagram is the coefficient with which the 
corresponding integral, defined according to the rules reviewed in fig. [IJ enters into our result for 
M^' D ~ 4 . An overall factor of 1/16 is suppressed and it is understood that one should sum over 
the 12 cyclic and reflection permutations of the external legs. In each coefficient, the second factor 
is a symmetry factor that accounts for overcounting in this sum. 
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Moreover, an integral appears with coefficient zero if and only if the integral is ill- 
defined (unregulated) after taking its external legs off shell and taking e — > 

Finally, it has been proposed that the signs ±1 of the contributing integrals can be 
understood by the requirement of cancelling unphysical singularities 



It is clear from fig. [7J that the first of these does not hold for our representation of the six- 
point amplitude; in particular, some relative weights are ±2, and there is one weight of —4. 
An examination of the integrals in fig. m shows that the second observation also requires 
some modification. This should not be too surprising since at six points we expect that some 
of the well-defined integrals appear in non-MHV amplitudes but not in the MHV ones. It 
would certainly be very interesting to determine whether any of these considerations could 
be generalized or modified to explain the pattern of coefficients appearing in fig. [7J 

IV. RESULTS 

Because the two-loop iteration formula (12.41) incorporates the known infrared singular- 
ities, it must hold for infrared-singular terms. We have evaluated the integrals in figs. [3] 
and H] through 0(e~ 2 ) analytically in terms of ordinary polylogarithms. In appendix lAl we 
have collected their values though 0(e~ 3 ); we refrain from presenting the much lengthier 
0(e~ 2 ) contributions. By inserting the values of the two-loop integrals into the assembly 
equation (13. 9p we find that eq. (12 .4p holds analytically through 0(e~ 2 ). This provides a non- 
trivial check on our cut construction, evaluation of integrals, and assembly of contributions. 
The reader may check agreement through (D(e~ 3 ) using the values of the two-loop integrals 
in appendix El and the one-loop amplitude (12. lOf) with n = 6. 

Beyond 0(e~ 2 ), we resort to numerical integration. We first constructed Mellin-Barnes 



representations, in order to make use of the package MB [27]. The package AMBRE 26] pro- 
vides a simple means for obtaining Mellin-Barnes representations that can be integrated 
using MB. One must treat the most complicated double-pentagon integrals, J*- 12 -* and I^ 13 ', 
with some care to produce a numerically suitable representation, so we give Mellin-Barnes 
representations for these two integrals in appendix O 

In four dimensions, there are at most four linearly independent momenta. For the six- 
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point amplitude, therefore, the Gram determinant of any five external momenta must vanish, 

det(Jfei • kj) = , i,j = 1,2,3,4,5. (4.1) 

This constraint turns out not to be relevant for any of our checks, but it is important to 
choose at least a few kinematic points satisfying this constraint, in order to ensure that 
any deviation from the BDS ansatz does not arise from choosing momentum invariants that 
cannot be realized in four dimensions. 



A. The /i-Dependent Terms 

In section II I II we split both the one- and two-loop amplitudes into a D = 4 part and a 
part containing explicit dependence on /z, the (— 2e)-dependent part of the loop momentum, 
according to eqs. (13.11) and (13.91) . The /i-dependent part of the one-loop amplitude, Mg M (e), 
vanishes as e — > 0. It only contributes to the remainder function R ( 6 , defined in eq. (12.91) . 
because it appears in (M^ 1} (e)) 2 multiplied by the singular terms in the one-loop amplitude, 
which are given in eq. (I2.10p . Thus the contribution of the one- and two- loop /z-dependent 
terms to is 



R 



(2),M 



lim 



1 1 

2^ 



i=i 



(4.2) 



where we set the dimensional regular izat ion scale \i — *■ 1 here to avoid confusion with \i v . 

Integrals containing numerator factors of fi p and \i q can be computed by differentiating 
integrands for scalar integrals with respect to Schwinger parameters 5?j]. This result holds 
because the dependence of the integrals on the (— 2e) components of the loop momenta is 
very simple. At two loops, it is given by, 



d fip CI fig exp [ip 1 p jjjq 1 g fJjp_^_q -L pq 



oc A e 



(4.3) 



where T p , T q and T pq are the sums of Schwinger parameters for propagators carrying loop 
momenta p, q and p + q, respectively, and A = T p T„ + T p T pq + T q T pq . Differentiation leads 
to the parameter insertions (see eq. (4.26) of ref. [5T ] ) . 
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f'p 



— 6- 



rp , rp 
- 1 q ~ ± pq 

A 



-e(T p + T q ) - (T q + T pq ) - (Tp + Ti 



m) 



T 
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(4.4) 
(4.5) 
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At one loop, jJl — > — e/T, where T is the sum of all the Schwinger parameters. 

At one loop, the insertion of 1/T shifts the dimension of the integral from D = 4 — 2e to 
D = 6 — 2e, which makes the integral infrared finite (and it remains ultraviolet finite). Thus 
J hex (e), and hence M^'^(e), vanish as e — > 0. At two loops, the factor of 1/ A in eqs. (I4.4p 
and (14.51) also shifts the dimension to D = 6 — 2e. However, the Schwinger parameters in the 
numerator lead to doubled propagators (see, for example, the discussion in ref. [57]]), which 
can cause infrared divergences, even near D = 6. 

In the case of the double pentagon integral eq. (14.51) shows that the doubled propa- 
gator is the central one. Because this propagator does not touch any on-shell external legs, 
doubling it is "safe", and we expect that i^ 14 -* = 0(e). We have confirmed this expectation, 
both numerically, and by checking analytically the analogous planar double box integral for 
massless four-point kinematics. 

In the case of the "hexabox" integral, I^ 15 \ eq. (14.41) leads to doubled propagators on the 
box loop, which create infrared divergences. However, the hexagon loop remains infrared 
safe. Reasonin g b y analogy to the factorization of soft and collinear singularities at the 



amplitude level 



281 ]. the /J?- hexagon inside the hexabox can be thought of as a hard process. 



Thus one can shrink it to zero size in space-time, and decorate it by a one-loop scalar triangle 
integral representing the infrared- divergent contributions. Thus we expect, 

I (15) (e) = -^(-sei)- 1 " 6 /^(e) + 0(e) . (4.6) 

The reason the equation is valid to 0(e) rather than O(e ) is simply because the "hard" 
part J hcx is itself 0(e). Once again, we have checked eq. (14.61) numerically. We have also 
checked analytically that the same relation holds for the analogous planar double box inte- 
gral. Inserting eq. (14.61) into the explicit expressions for the one- and two-loop //-dependent 
contributions Mg 1 ^ and M^'^ in eq. (14.21) . we see that R^'^ vanishes. 



B. Evaluation of Remainder Function 



Explicit computations 



14 



15] have demonstrated that the remainder function R, 



(2) 



defined in eq. (12. 9ft vanishes for n — 4, 5. In this section we shall evaluate R 6 numerically 
at a few kinematic points, and find that it is nonzero and nonconstant. 

We choose Euclidean kinematics for all points, as this simplifies the numerical evaluation. 
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A particularly convenient kinematic point is 

: Si t i+\ = —1 , Si,i+i,i+2 = —2 , 



(4.7) 



which we take to be our standard reference point. This point has several advantages. Firstly, 
because it is symmetric under cyclic relabeling % — > i+1 and under the reflection i — > 6— i+1, 
we do not need to evaluate any relabelings of the integrals to obtain the amplitude. We 
have also exploited s^i+i = — 1 to simplify the Mellin-Barnes representations. Moreover, 
this kinematic point satisfies the Gram determinant constraint (14.11) . 

The numerical values of all the integrals in fig. [3] at the standard kinematic point are given 
in appendix [B] (For completeness we also give the values of the non-contributing integrals 
in fig. HI) Inserting these values into the amplitude (I3.10p . we obtain, 

Mf < D = 4 = A - - ^ - 20.8534 ± 0.0057 + 0(e) , (4.8) 

at the standard kinematic point K^°\ The values of the /i-dependent contributions (13. lip 

are, 

M 6 (2) ' M = + 8.6024 ± 0.0010 + 0(e) . (4.9) 



We have included estimated errors from the numerical integration reported by CUBA [59], 
added in quadrature. In many cases the errors appear to be overestimated. However, 
there can be correlations between different subintegrals in which the amplitude is expanded, 
because the same random seed is used. The reported errors do appear to give us a reliable 
measure of how many digits are trustworthy. 

The values (14. 8 j) and (14.91) may be compared to those for the ABDK/BDS ansatz (12.41) . 
Separating it out, in analogy to ( 13. 91) . as 

M 6 BDS = M 6 BDS ' D=4 + M 6 BDS '" , (4.10) 

where the second term arises from J hex in Mg , we find that 

= £ - ^ - ^ _ 21 . 947 1 + 0(e) , (4.11) 



2e 4 e 
and that the \i pieces are given by 



m bds,^ = 2^3510 + g goi7 + 0(e) (4 J2) 
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Since these formulae involve computing only one-loop integrals, the numerical integration 
errors are much smaller and do not affect the answer to the quoted precision. 

By comparing eqs. H4.9[) and (I4.12p . we see that the /i terms agree. This result is in accord 
with the general vanishing of R^'^ described in section HV Al However, there is a difference 
in the D = 4 terms, between our explicit calculation of the amplitude and the BDS ansatz. 
Defining it as 

R A = Rf = M 6 (2) - M 6 BDS , (4.13) 
we find that at our standard kinematic point (14.71) it equals 

R° A = R A (K {0) ) = 1.0937 ± 0.0057 . (4.14) 

Although the remainder is only 5 percent of the finite term in eq. (14.81) . it is nonzero at very 
high confidence level, demonstrating that the ABDK/BDS ansatz needs to be modified. 

Besides our standard kinematic point, we also evaluated Ra at various other kinematic 
points, 



K (l) 


■ «12 = 


-0.7236200, 


s 23 = -0.9213500, 


S34 = - 


-0.2723200, 


S45 = 


-0.3582300, 




«56 = 


-0.4235500, 


s 6 i = -0.3218573, 


S123 = 


-2.1486192, 


S234 


= -0.7264904, 




S345 = 


= -0.4825841, 














: S12 = 


-0.3223100, 


s 23 = -0.2323220, 


S34 = - 


-0.5238300, 


S45 = 


-0.8237640, 




s 56 = 


-0.5323200, 


s 61 = -0.9237600, 


S123 — 


-0.7322000, 


•^234 


= -0.8286700, 




S345 - 


= -0.6626116, 














: Sj^+i 


= -1, S123 = 


-1/2, S234 = -h/i 


s, S345 


= -17/14, 








: Sj^+i 


— — 1, Si,i+l,i+2 — — 3, 












: s i,i+l 


= —1? Si,i+l,i+2 — —9/2. 








(4.15) 



With six-point kinematics we have sufficient freedom to construct three nontrivial con- 
formal cross ratios: 



x 13 x 46 


S12S45 




S123S345 


,7^2 ™2 


S23S56 


/-y»2 ,-y»2 

• L 25 Jj Al 


■5234-5123 


,7^2 ,y,2 


S 34 S 6 i 


/yi2 ,71 2 

x 36 x 52 


•S345S234 



(4.16) 
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TABLE I: The numerical remainder compared with the ABDK ansatz (12. 4h for various kinematic 
points. The second column gives the conformal cross ratios defined in eq. (|4.16p . 



kinematic point 


(ni,n 2 ,n 3 ) 


Ra 


K (0) 


(1/4,1/4,1/4) 


1.0937 ±0.0057 




(1/4,1/4,1/4) 


1.076 ± 0.022 




(0.547253, 0.203822, 0.881270) 


-1.659 ±0.014 




(28/17,16/5,112/85) 


-3.6508 ± 0.0032 


K w 


(1/9,1/9,1/9) 


5.21 ±0.10 


K(5) 


(4/81,4/81,4/81) 


11.09 ±0.50 



The three conformal cross ratios (14.161) for these kinematic points are given in the second 
column of Table [B Here K^ l > , and K^ 3 ' satisfy the Gram determinant constraint (14. ip 
while if (4) and K® do not. 

The point K^ 1 ' is chosen so that the conformal cross ratios Ui = 1/4 are identical to 
the cross ratios for our reference kinematic point K^°\ The agreement, within the errors, 
between the remainder functions for these two kinematic points, suggests that Ra is a 
function of only the cross ratios, i.e. is invariant under dual conformal transformations. 



C. Comparison with Wilson Loop 



Drummond, Henn, Korchemsky and Sokatchev have already shown 20J that the Wilson 

(2) 

loop expectation value (Wg ) corresponding to the two- loop six-point MHV amplitude, is 
not equal to that suggested by the ABDK/BDS ansatz for amplitudes. That is, they found 
a nonvanishing remainder function, 



(4.17) 



Here W^ DS is the Wilson loop analog of Mq DS defined in eq. (12. 6p . It involves the same 
function Mg^e), but different constants appear {20], [21 j than the ones for amplitudes (which 
are captured by /^(e) and C®). We are motivated by the correspondence between MHV 
amplitudes and Wilson loops at one loop (si, 9], to ask how these two remainder functions, 
Ra and Rw, compare. 
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TABLE II: The comparison between the remainder functions Ra and R\y for the MHV amplitude 
and the Wilson loop. To account for various constants of the kinematics, we subtract from the 
remainders their values at the standard kinematic point K(°\ denoted by R A and R w - The third 
column contains the difference of remainders for the amplitude, while the fourth column has the 
corresponding difference for the Wilson loop. The numerical agreement between the third and 
fourth columns provides strong evidence that the finite remainder for the Wilson loop is identical 
to that for the MHV amplitude. 



kinematic point 




Ra — R\ 


Rw — R%/ 




(1/4,1/4,1/4) 


-0.018 ±0.023 


< 10" 5 


K m 


(0.547253, 0.203822, 0.881270) 


-2.753 ±0.015 


-2.7553 


K m 


(28/17, 16/5,112/85) 


-4.7445 ± 0.0075 


-4.7446 


#(4) 


(1/9,1/9,1/9) 


4.12 ±0.10 


4.0914 




(4/81,4/81,4/81) 


10.00 ± 0.50 


9.7255 



As explained in ref. [20(, in the collinear limits, corresponding to u\ = 0,us = 1 — u 2 , 
the Wilson loop remainder function Rw becomes a constant, corresponding to eq. (17) that 
paper. On the other hand, as explained in section |V] the MHV amplitude remainder function 
must vanish in the collinear limits in order to be consistent with collinear factorization. This 
suggests a simple relation between the two remainders, 



Ra — Rw — cw , 



which we test numerically. From DHKS 21), the constant cw takes on the value, 

c w = 12.1756, 



(4.18) 



(4.19) 



with a precision of ~ 1CT 3 . 

The numerical determination of the constant cw from the collinear limits of the Wilson 
loop leads to some loss of precision, so instead in Table [TT] we compare the Wilson loop 
and MHV-amplitude remainder functions by considering a "difference of differences" . That 
is, for both the Wilson loop and MHV amplitude the remainders Ra and Rw are found 
by subtracting the value of the appropriate ABDK/BDS formula at that point. From Ra 
and Rw we subtract the corresponding values R° A and R w at the standard kinematic point 



31 



From DHKS the remainder at the standard point is 2l| = 13.26530. This 
subtraction eliminates any dependence on cw- The Wilson loop results are obtained from 



ref. 



2l| . In general, the Wilson loop results have much smaller errors than those of the 



amplitude. This is due to the much simpler integral representations appearing in the Wilson 
loop computation [211 ]. 

Using the value of cw from eq. (j4.19j) . we can also compare Ra and Rw directly, albeit at 
lower precision. We find that within errors eq. (14.181) is satisfied for all six kinematic points 
in eqs. (14. 7ft and ( 14.151) . As mentioned above, we tested dual conformal invariance directly 
at one point. However, this invariance was tested much more extensively (also numerically) 



for the Wilson loop 



20 



211 ] . Thus our numerical agreement with the Wilson loop remainder, 



displayed in Table [TTl obviously provides considerable additional evidence that Ra possesses 
dual conformal invariance. 



V. THE REMAINDER FUNCTION 

In the previous section we found a numerical difference between the ABDK/BDS ansatz 
for the two-loop six point amplitude and the explicit calculation. In this section we constrain 
the analytic form of the remainder and point out that its functional form can be determined 
from triple-collinear limits. 

A. Constraints on the Remainder Function 

The form of the ABDK/BDS ansatz is tightly constrained by factorization properties 
and also exhibits dual conformal invariance. As discussed in the previous section, numerical 
evidence confirms that, while it departs from the ansatz, the even part of the two-loop six- 
point amplitude is invariant under dual conformal transformations. We can therefore discuss 
further constraints imposed by this symmetry on the remainder function. We will discuss 
the case of n-particle amplitudes and specialize to n = 6 at the end. 

DHKS [3] argue that MHV amplitudes (like Wilson loops), should obey anomalous dual 
conformal Ward identities, the anomaly due to infrared (ultraviolet) divergences. The fact 
that the BDS ansatz accounts for all infrared divergences of MHV amplitudes to all loop 
orders implies that R n is finite and thus independent of the regulator. It also means that R n 
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should satisfy non-anomalous Ward identities; that is, it must actually be invariant under 
dual conformal transformations and thus depend only on conformally-invariant cross ratios: 



Rn = Rn({Uijkl}) • (5.1) 

Here (ijkl) denote the allowed quartets of the external legs leading to well-defined cross ratios 
(12.3ip . The symmetry properties of the remainder function Rn under the permutation of 
its arguments follow from the reflection and cyclic identities obeyed by the rescaled MHV 
scattering amplitudes M n L \ 

Further restrictions on R n come from the fact that the BDS ansatz correctly captures 
the two-particle collinear factorization of MHV amplitudes. In general, the L-loop rescaled 



planar amplitudes M n L \l, 2, . . . , n) satisfy simple re 



adjacent legs k iy k i+1 become collinear, |32|, |60|, |6l|, |62| 



ations as the momenta of two color- 



Mf )(. . . , i*% (,• + i) A ^, . . o — E E ri ^ {z] iXi > {t + 1 ) Ai+1 ) M -i°(- • • * p " • • •) • ( 5 - 2 ) 

1=0 x=± 

The index I sums over the different loop orders of the rescaled splitting amplitudes, 



r%{z; i\ (i + 1)^) = 7"7"r' 1 ; , (5.3) 



Split W A (z;z A ',(i + l) A ^; 
SplitL°i(^;i A s(i + l) Ai+i ; 



while A sums over the helicities of the intermediate leg kp = (hi + k i+ i), and z is the 
longitudinal momentum fraction of fcj, fcj w zkp. 

The relevant two-loop splitting amplitudes were calculated in refs. (3), [s^]. If we assume 
that dual conformal symmetry holds to all orders, then the five-point amplitudes are fully 
determined by this symmetry. By taking the collinear limit, the all- loop splitting amplitude, 



r fuii 



l + ^a L r (i) , (5.4) 

L=l 

must have the form 

00 

In r fu11 = E a ' / (0 ( £ ) r(1) ( /e ) + 0{e). (5.5) 

n 

As discussed in ref. |6[, this iterative structure yields the correct collinear behavior to all 
loop orders. 

Since the BDS ansatz accounts for collinear factorization, to all orders in perturbation 
theory, the remainder functions must have a trivial behavior under collinear factorization. 
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The n- and (n — l)-point remainder functions must be related by 

lim Rn({u iu i 2 ,i 3M }) = Rn-l({Ui u i 2 ,i 3M Y) , (5.6) 



for any two-particle Mandelstam invariant xj i+2 = Sj,i+i- The arguments of R n -\ are the 
subset of the n-point conformal cross ratios that are non-vanishing and well-defined in the 
collinear limit s^i+i — > 0. 

For five-point amplitudes, no conformal cross ratio with the required properties may be 
constructed. Thus no remainder function can exist consistent with collinear factorization 
and the requirement of dual conformal invariance. A constant remainder is ruled out by 
collinear factorization. 

Nontrivial conformal cross ratios can be first constructed with six-particle kinematics, 
«i, u 2 and u 3 in eq. (14.161) . The cyclic and reflection symmetries imply that R e (ui,u 2 ,u 3 ) is 

(2) 

a totally symmetric function of its arguments. Because the remainder function R 5 for the 
five-point amplitude vanishes Q, [3], R^ (ui,U2,Us) must vanish in all collinear limits. If 
dual conformal symmetry is valid to all loop orders then i?5 vanishes exactly and therefore R$ 
must vanish in all collinear limits to all loop orders. DHKS [^(J reached a similar conclusion 
in their analysis of the two-loop six-sided Wilson loop, namely that there is a remainder / 
above and beyond the ABDK/BDS ansatz with the properties described above. 



B. Remainder Function from Triple- Collinear Limits 

Planar color-ordered scattering amplitudes exhibit singularities when several adjacent 
momenta become collinear. General all-order factorization properties of scattering ampli- 
tudes have been discussed in ref. j^. The most familiar of these limits are when just two 
particles become collinear. However, the limits when more particles become collinear simul- 
taneously can provide additional constraints 4 . (Multi-collinear configurations should not be 
confused with multi-particle factorization limits, in which amplitudes factorize into products 
of lower-point, non-degenerate scattering amplitudes. The latter limits are trivial for MHV 
amplitudes in supersymmetric theories.) 

As discussed previously, the ABDK/BDS ansatz incorporates the correct two-particle 
collinear factorization of MHV amplitudes. It also makes definite predictions, which remain 

4 We thank Gregory Korchemsky and Emery Sokatchev for discussion on the triple-collincar limits. 
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to be tested, for the multi-collinear factorization of MHV amplitudes in M = 4 super- Yang- 
Mills theory. These tests amount to constraints on the remainder functions R n . For the 
six-gluon amplitude we only have the triple-collinear limit. As we will see, it is possible to 
completely determine the remainder function P 6 from this limit. 

Let us consider three adjacent momenta fc a ,6,c hi the limit that they become collinear and 
introduce the three momentum fractions 

k a = z\P , k b = z 2 P , k c = z 3 P , zi + z 2 + z 3 = 1 , < Zi < 1 , P 2 -> . 

(5.7) 

An n-point amplitude at I loops factorizes as follows: 

l 

- ^2J2 A n-2\h,...,k n - 3 ,P x ) S V lit ( ll(k n ^k n ^k n ;P). 

A=± s=0 

(5.8) 

Taking into account parity and reflection symmetries, there are six independent triple- 
collinear splitting amplitudes: 

Sj>lit + {k+k*k+;P), (5.9) 
Split_ Ap (A;^A; c Ac ;P), A a + \ b + A c - A P = 2 , (5.10) 

Split_ Ap (^A; fe Hh P), A a + X b + A c - X P = . (5.11) 

The first one (15.91) vanishes in any supersymmetric theory. The three triple-collinear splitting 
amplitudes of the second type (15.101) . an example of which is A a = A& = A c = Xp — 1, 
appear in limits of MHV amplitudes. The Af = 4 supersymmetry Ward identities for MHV 
amplitudes imply that their rescaled forms 5 are all equal, 

SpUtg^fc+fc+jP) Split^(fc+fc^+;P) (t) Shr 
S P litf(ktk+kf,P) SplitP(k+k±ki;P) S ^^ U3) - 
The two splitting amplitudes of the third kind (15. lip arise only in limits of non-MHV 
amplitudes and do not have a simple factorized form similar to (15. 121) . 6 



5 We omit a trivial dimensional dependence on s a & c from the argument list of rl 



(> 



The spin-averaged absolute values squared of tree-level triple-collinear splitting amplitudes have been 



computed in ref. 6J]; without spin-averaging they have been computed in refs. 65(. The tree- level triple 



(and higher) collinear splitting amplitudes themselves have been computed in ref. 66|, and in ref. 67 1 



using the MHV rules 68|. The one-loop correction to the q — * qQQ triple-collinear splitting amplitude in 



QCD was computed in ref. 69]. 
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On general grounds, the six-gluon amplitude exhibits a nontrivial triple-collinear limit. 
In the limit (|5.7p . for a = 4,6 = 5,c = 6, the three conformal cross ratios (14.161) are all 
nonvanishing and arbitrary: 

■545 1 _ S 56 1 _ Z\Z 3 

Ul = , u 2 = , U 3 = — r. (5.13) 

5456 1 -^3 «456 I ~ Z 1 (1 - Z X ) (1 - Z 3 ) 

The remainder function therefore survives the triple-collinear limit, and is evaluated at 
(ui, U2, u 3 )- Thus, assuming dual conformal invariance, finding the remainder function in 
this limit is equivalent to finding it for generic six-particle kinematics. 

We can expose the two-loop remainder function by comparing eq. (15. 8ft for the rescaled 
MHV six-point amplitude, 

lim M 6 (2) = Mf ) +M] 1) r ( t ! ) (^,^,z 1 ,Z3,e)+rl ; 2) (^,^,z 1 ,Z3,e), (5.14) 

4 || 5 ||g b 4 1 4 6 V S456 ' S456 ' L ' a ' 1 b V S456 's 45 6 ' 

with the triple-collinear limit of eq. (12.91) . The second term in eq. (15.141) is already determined 
by one-loop calculations, and it is incorporated in the ABDK/BDS ansatz. Therefore -Rg 2 '* 
enters only in the two-loop splitting amplitude rg , as the deviation from the ABDK/BDS 
prediction r^ BDS : 

r s ] (^ *> e ) = r ? )BDS (^' e ) + fia, «3) , (5-15) 



with 



r (2)BDS^ ^ x 



1 ( m, ... ... \ 2 

2 

Thus, the two-loop remainder function is completely determined by the two-loop triple- 
collinear splitting amplitude, e.g for the helicity configuration (15.101) . While a priori it may 
depend on all four arguments of the splitting amplitude, dual conformal invariance requires 
that it depend only on the three cross ratios. 

The triple-collinear splitting amplitudes may be computed using the unitarity method 



following the strategy in refs. 62|, |63J]. It is important to understand whether they satisfy an 
iteration relation generalizing that of the double-collinear splitting amplitude 3] . If such an 
iteration relation exists, then it should be straightforward to construct an all-order iteration 
relation for the six-point gluon amplitude. This would allow us to add in a correction term 
to the BDS ansatz, at least for the six-point case. 
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The remainder function beyond two loops can also be extracted from the triple-collinear 
splitting amplitude, though they are no longer equal. Instead, they are related iteratively 
via, 

i 

s=2 

For amplitudes with additional external legs, it is unclear whether the triple-collinear 
limits suffice to constrain the remainder functions completely. If these limits do not suffice, 
we can formulate additional constraints along the lines above. In particular, it is easy 
to see that the remainder function of the two-loop n-point MHV amplitude is completely 
determined by the difference between the two-loop (n — 3)-point splitting amplitude and 
the iteration of the one-loop {n — 3)-point splitting amplitude. Also, the consideration of 
the m-particle collinear limit with m < n — 4 leads to consistency conditions analogous to 
eq. (ED. 



■^S in- zuz 3 ,e) 

S456 ' S 45 6 ' X ' J ' ' 



(s)BDS, S45 
' ' S 1 



\ 5456 ' «456 ' ' 3 ' ' 



. (5.17) 



VI. SUMMARY AND CONCLUSIONS 

The M = 4 supersymmetric gauge theory has proven an important laboratory and testing 
ground for inquiry into the properties of gauge theories, both at weak and strong coupling. 
The BDS ansatz for planar MHV scattering amplitudes p] in this theory, along with the BES 
integral equation for the cusp anomalous dimension [3] and the strong-coupling calculation 
of Alday and Maldacena pj] , point to the possibility of computing planar amplitudes for any 
value of the coupling. 

In this paper we have checked the BDS ansatz directly by computing the parity-even 
parts of the leading-color part of the planar two-loop six-point amplitude. Using the uni- 
tarity method, we have obtained an integral representation for it. Numerical evaluation of 
this representation shows that there is a remainder beyond the ABDK/BDS prediction for 
this amplitude. Strikingly, the remainder agrees with the corresponding remainder for the 



hexagonal Wilson loop 2Cj, [21 ] . 



This remainder must vanish in any limit where two color-adjacent momenta become 
collinear, because the ABDK/BDS construction accounts for all terms with collinear singu- 
larities. As we showed, it should be possible to fully reconstruct the remainder function for 
the six-point amplitude by evaluating triple-collinear splitting amplitudes. 
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There are a number of interesting open issues which remain to be clarified. The origin 
of the dual conformal symmetry remains mysterious. In the context of the AdS / CFT corre- 
spondence it has been suggested [16| that it is related to symmetries of the space defined by 
the coordinates introduced in eq. (12.24)) . Based on this, one might wonder whether dual 
conformal symmetry can be found in the planar amplitudes of all four-dimensional CFTs 
with a string-theory dual. 

We have found that the integrals appearing in the amplitude either vanish as the loop mo- 
menta are taken to be four dimensional or are pseudo-conformal (with the external momenta 
taken off shell to make them infrared finite). Contributions containing the (— 2e)-dimensional 
components of loop momenta in the numerator factors satisfy the BDS ansatz and drop out 
of the remainder function. It seems reasonable to expect that this pattern continues to all 
loop orders. It would therefore be very useful to have a set of rules for writing the coefficients 
of all pseudo-conformal integrals directly, without resorting to evaluations of the cuts. 

So far dual conformal transformations have only been discussed in the context of their 
action on Wilson loops or planar MHV amplitudes after dividing out the tree amplitude 
prefactor. Can we extend this to non-MHV amplitudes? At least at one loop the integrals 
appearing in non-MHV amplitudes are pseudo-conformal scalar box integrals 32j, hinting 
that dual conformal invariance might be a general property of the planar limit of the theory. 
However, in this case the tree amplitude does not factor out [45j, |46] , leaving the question of 
how the dual conformal symmetry might act on the spinor products that enter the relative 
factors of different pseudo-conformal integrals. Related to this is the question of whether 
the dual conformal symmetry can be extended to the Lagrangian of M = 4 super- Yang-Mills 
theory. It is not clear how it should act on the Lagrangian, given that it is only understood 
at present for planar MHV amplitudes. 



What other properties can we employ to constrain the scattering amplitudes? Integra- 



bility 
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70j of the dilatation operator in planar M = 4 super- Yang-Mills theory has 



not yet been used. Similarly, we expect the amplitudes to have simple structures in twistor 
space 711 ] . At one loop the coefficients in front of each integral have been shown to lie on 
simple curves in twistor space, though the precise structure of complete loop amplitudes has 



not been determined 
higher loops. 



46, 



53, 



721 ] . It would be very interesting to explore the structure at 



In order to shed light on the remainder function it is very important to find its analytic 
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form for the two-loop six-point amplitude. This form would be extremely useful for under- 
standing the missing terms in the BDS ansatz at higher loops. It would also be useful for 
analytic continuation into the physical high energy or Regge limits of 2 — > 4 and 3^3 scat- 



tering that have been discussed recently 



18 



19| . In section |V] we proposed a possible means 



for constructing it at six points from triple-collinear limits. The equality of the six-point 



remainder function and the corresponding Wilson loop quantity [2JJ is rather surprising. 
It would obviously be very desirable to understand whether this equality holds to all loop 
orders, as well as at strong coupling. 

In summary, our computation demonstrates that the BDS ansatz requires modification 
for amplitudes with six or more external legs. The surprising equality of the Wilson loop 
and MHV amplitude remainders, however, points to an additional structure in the theory 
which constrains its form. This in turn provides hope of determining the remainder function 
analytically, first at two loops, and eventually to all loop orders. 
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APPENDIX A: ANALYTIC VALUES OF INTEGRALS 



1. One-loop integrals 



In the Euclidean region with all s^j+i and s iji+ i )i+ 2 negative, the one-mass box, I lm in 
fig. El is given through O(e ) by 
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Similarly, the two-mass "easy" box, 7 , is 
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2. Two- loop integrals 



We have computed analytic expressions for the two-loop integrals appearing in fig. [3] 
through 0(e~ 2 ). The 0(e~ 2 ) expressions are rather cumbersome, so we display here only 
the results through (9(e -3 ), omitting as it is given by a product of two one-loop integrals 
of type J lm : 
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APPENDIX B: NUMERICAL VALUE OF INTEGRALS 



In this appendix we give numerical values of the one- and two-loop integrals at the 
symmetric kinematic point K^°\ 

The numerical values of the one-loop integrals flAlj) and (1A2|) are needed through 0(e 2 ), 
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The one-loop hexagon in fig. [2] is 
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If no errors are quoted, the integration errors are smaller than the quoted precision. 



APPENDIX C: MELLIN-BARNES REPRESENTATIONS 



(B3) 



In this appendix we present Mellin-Barnes representations of the most complicated inte- 
grals of fig. El 7 (12) and J (13) . In both cases it is necessary to introduce an auxiliary parameter 
7] in order to render the integral well-defined. The value of the integral in the desired limit 
7] — > is obtained by analytic continuation. Using the notation Zi j = Zi + Zj + • • ■ we have: 
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